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1. Introduction 

With the recent announcement [1] of an Englert-Brout-Higgs (EBH) [2] candidate boson af- 
ter the start-up and successful running of the LHC, the era of precision QCD, wherein the total 
precision tag is 1% or better, is upon us. The attendant need for exact, amplitude-based resumma- 
tion of large higher order effects is now more paramount, in view of the expected role of precision 
comparison between theory and experiment in determining the detailed properties of the newly dis- 
covered EBH boson candidate. It has been argued elsewhere [3,4] that such resummation allows 
one to have better than 1% theoretical precision as a realistic goal in such comparisons. Here, we 
present the status of this approach to precision QCD for the LHC with the attendant IR-improved 
DGLAP-CS [5,6] theory [7, 8] realization via HERWIRI1.031 [9] in the HERWIG6.5 [10] envi- 
ronment in interplay with NLO exact, matrix element matched parton shower MC precision issues. 
We employ the MC@NLO [11] methodology to realize the attendant exact, NLO matrix element 
matched parton shower MC realizations in comparisons with recent LHC data. 

In the discussion we continue the strategy of building on existing platforms to develop and 
realize a path toward precision QCD for the physics of the LHC. We exhibit a union of the new 
IR-improved DGLAP-CS theory and MC@NLO. We are also pursuing the implementation [12] of 
the new IR-improved DGLAP-CS theory for HERWIG++ [13], HERWIRI++, for PYTHIA8 [14] 
and for SHERPA [15], as well as the corresponding NLO ME/parton shower matching realizations 
in the POWHEG [16] framework - one of the strongest cross checks on theoretical precision is the 
difference between two independent realizations of the attendant theoretical calculation. 

We set the stage for the proper exposition of the interplay between the NLO ME matched 
parton shower MC precision and the new IR-improved DGLAP-CS theory in the next section by 
showing how the latter theory follows naturally in the effort to obtain a provable precision from 
our approach [4] to precision LHC physics. We review this latter approach in the next section as 
well. We then turn in Section 3 to the applications to the recent data on single heavy gauge boson 
production at the LHC with an eye on the analyses in Refs. [9] of the analogous processes at the 
Tevatron. We will focus in this discussion on the single Z/y* production and decay to lepton pairs 
for definiteness. The other heavy gauge boson processes will be taken up elsewhere [12]. 



2. Brief Recapitulation 

The starting point for our discussion may be taken as the fully differential representation 

do = ^ / dx\dx2Fj(xi)Fj(x2)dd res (xiX2s) (2.1) 

of a hard LHC scattering process using a standard notation so that the {Fj} and <ia res are the re- 
spective parton densities(PDFs) and reduced resummed hard differential cross section, where the 
resummation is for all large EW and QCD higher order corrections in order to achieve a total pre- 
cision tag of 1% or better for the total theoretical precision of (jO|). The proof of the correctness of 
the value of the total theoretical precision Aa t h of Q2.1| ) is the basic issue for precision QCD for the 
LHC . This precision can be represented as follows: Aa t h = AF©Aa res where AA is the contribution 
of the uncertainty on the quantity A to AoVh 1 - The proof of the correctness of the value of the total 
theoretical precision Aa t h is essential for validation of the application of a given theoretical predic- 
tion to precision experimental observations for the signals and the backgrounds for both Standard 
Model(SM) and new physics (NP) studies, and more specifically for the overall normalization of 
the cross sections in such studies. We cannot emphasize too much that NP can be missed if a calcu- 
lation with an unknown value of Aa t h is used for the attendant studies. We note that here Aa t h is the 
total theoretical uncertainty that comes from the physical precision contribution and the technical 

precision contribution [17]: the physical precision contribution, Aa£ hys , arises from such sources 
as missing graphs, approximations to graphs, truncations,....; the technical precision contribution, 
Aa t ^ ch , arises from such sources as bugs in codes, numerical rounding errors, convergence issues, 
etc. The total theoretical error follows from 

Ac7 th = Aa t p h hys ©Aar h . (2.2) 



Here, we discuss the situation in which the two errors in the equation for Aa^ are independent for definiteness; the 
equation for it has to be modified accordingly when they are not. 
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As a general rule, the desired value for Aa t h, which depends on the specific requirements of the 
observations, should fulfill Aa t h < /Aa expt . Here Aa expt is the respective experimental error and 

/ < \ so that the theoretical uncertainty does not significantly adversely affect the analysis of the 
data for physics studies. 

With the goal of realizing such precision in a provable way, we hav e developed the QCD 



QED resummation theory in Refs. [4] for the reduced cross section in ( 24 ) and for the resummation 
of the evolution of the parton densities therein as well. Mainly because the theory in Refs. [4] is 
not widely known, we recapitulate it here briefly. The master formula for our resummation theory 
may be identified as 



_ „SUMir(QCED) v ~ _J_ f U n d%_ 
""res — e Ln,m=0 n\m\ J 11/1=1 k h 

nm d3k 'i 2 f _d*y Jy-(pi+q[-p2-q2-Y,kj -Y.^ h)+DQ0SD 
72 = 1 k' h J (27rf e 

fin,m{k\ j • • • i^n'i^\ > • ■ • j k m ) J? i (2-3) 



where <ia res is either the reduced cross section dd res or the differential rate associated to a DGLAP- 
CS [5,6] kernel involved in the evolution of the {Fj} and where the new (YFS-style [18]) non- 

Abelian residuals fi nm (k\ , . . . ,k n ;kf l , . . . ,k' m ) have n hard gluons and m hard photons and we show 
the final state with two hard final partons with momenta p2, q2 specified for a generic 2f final state 
for definiteness. The infrared functions SUMtr(QCED), Dqced are defined in Refs. [4,7,8]. This 
simultaneous resummation of QED and QCD large IR effects is exact. 

We have shown in Refs. [7-9] that the methods in Refs. [19,20] give approximations to our 



hard gluon residuals j3„; for, the methods in Refs. [19,20], unlike the master formula in (2.3), are 
not exact results. The threshold-resummation methods in Refs. [19], using the result that, tor any 
function f(z), 



dzz n ~ l f(z) 



< (-) max \f{z) 



drop non-singular contributions to the cross section at z — > 1 in resumming the logs in «-Mellin 
space. The SCET theory in Refs. [20] drops terms of at the level of the amplitude, where X = 
^/K/Q for a process with the hard scale Q with A ~ .3GeV so that, for Q ~ lOOGeV, X = 5.5%. 
The known equivalence of the two approaches implies that the errors in the threshold resummation 
must be similar. We can only use these approaches as a guide to our new non-Abelian residuals as 
we develop results for the sub-1% precision regime. 

As it is explained in Refs. [4], the new non-Abelian residuals /3 m „ allow rigorous shower/ME 

matching via their shower subtracted analogs: /3„ u , — > ft mM where the p m n have had all effects in 
the showers associated to the {Fj} removed from them and this naturally brings us to the attendant 
evolution of the {Fj}. For a strict control on the theoretical precision in (p]), we need both the 
resummation of the reduced cross section and that of the latter evolution. 



When the QCD restriction of the formula in fl2.3[ ) is applied to the calculation of the kernels, 
Pab, hi the DGLAP-CS theory itself, we get an improvement of the IR limit of these kernels, an 
IR-improved DGLAP-CS theory [7-9] with new resummed kernels P^, which are reproduced 
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here for completeness: 

Kf{z) = C F F Y¥S {y q )e^ 



l-z 



(1 



f q (r q )8(i-z) 



P^(z)=C F F YFS ( Yq )e^ 1 + il Z) \ r«, 



2C g F yfs (7g)^ 5g { 



1 -z 



7g . 



l-z 



7c; 



+ 2(z 1+7G (l -z) + z(l -z) 1+ *=) -/ G (y G )5(l -z)}, 
P^Jiz) = F YFS (YG)e^{z 2 (l -z)* + (1 -z)V G }, 



(2.4) 



where the superscript "exp" indicates that the kernel has been resummed as predicted by Eq. ( [Q| ) 
when it is restricted to QCD alone, where the YFS [18] infrared factor is given by Fyfs(«T = 
e~ CE "/T(l +a) where Ce is Euler's constant and where the respective resummation functions 

Ya-, §a, Sa-, -A = q, G are given in Refs. [7, 8] 2 . Cp(Cc) is the quadratic Casimir invariant for the 
quark(gluon) color representation respectively. From these new kernels we get a new resummed 
scheme for the PDFs and the reduced cross section: 



Fj, &-+F}, a' for 
P Gq (z) 



Gq 1 



z ,etc, 



(2.5) 



with the same value for o in (|2JJ) with improved MC stability [9] - we do not need an IR cut-off 
'jfeo' parameter in the attendant parton shower MC based on the new kernels. Note that, while the 
degrees of freedom below the IR cut-offs in the usual showers are dropped in those showers, in 
the showers in HERWIRI 1.031, as one can see from (jO]), these degrees of freedom are integrated 
over and included in the calculation in the process of generating the Gribov-Lipatov exponents y& 
in (|2~!4|). The new kernels agree with the usual kernels at 0(oc s ) as the differences between them 
start in ^(q s 2 ), so that the NLO matching formulas in the MC@NLO and POWHEG frameworks 
apply directly to the new kernels for exact NLO ME/shower matching. 

In Fig. 1 we show the basic physical idea of Bloch and Nordsieck [23] underlying the new 
kernels: an accelerated charge generates a coherent state of very soft massless quanta of the re- 



'/ 



Soft Gluon Cloud 

E-E 



9(1 " ~) 



Gi, i = 0, ■ ■ ■ ,oo 



Figure 1: Bloch-Nordsieck soft quanta for an accelerated charge. 

spective gauge field so that one cannot know which of the infinity of possible states one has made 
in the splitting process q(l) — > q(l — z) + G(g> G\ ■■ ■ ® Ge, £ = 0,- ■■ ,°° illustrated in Fig. 1. The 

2 The improvement in Eq. ( b_.4| ) should be distinguished from the resummation in parton density evolution for the 
"z — ¥ 0" Regge regime - see for example Ref. [21, 22]. This latter improvement must also be taken into account for 
precision LHC predictions. 
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new kernels take this effect into account by resumming the terms G ( ( G£ s rn(^)ln(l — z)) n ) when 

z — > 1 is the IR limit. From (2.5) and (|Q|), one can see that when the usual kernels are used these 
terms are generated order-by^oraer in the solution for the cross section a in (2.1). Our resumming 
of these terms enhances the convergence of the representation in (2.1) for a given order of exactness 
in the in input perturbative components therein. In the next section we illustrate this last remark in 
the context of the comparison of NLO parton shower/matrix element matched predictions to recent 
LHC data. 



3. Interplay of NLO Shower/ME Precision and IR-Improved DGLAP-CS Theory 

Here, we compare new MC HERWIRI1.031 [9] with HERWIG6.510, both with and with- 
out the MC@NLO [11] exact ff(a s ) correction to illustrate the interplay between the attendant 
precision in NLO ME matched parton shower MC's and the new IR-improvement for the kernels 
realized in Herwiri 1.031, where we use the new LHC data for our baseline for the comparison. 

For the single Z/y* production at the LHC, we show in Fig. || in panel (a) the comparison 
between the MC predictions and the CMS rapidity data [25] and in panel (b) the analogous com- 
parison with the ATLAS Pj data, where the rapidity data are the combined e + e~ — }JL~ }JL + results 
and the pj data are those for the bare e + e~ case, as these are the data that correspond to the theoret- 
ical framework of our simulations - we do not as yet have complete realization of all the corrections 
involved in the other ATLAS data in Ref. [26]. These results should be viewed with an eye on our 
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Figure 2: Comparison with LHC data: (a), CMS rapidity data on (Z/y*) production to e + e , jU+jii 
pairs, the circular dots are the data, the green(blue) lines are HERWIG6.510(HERWIRI1.031); (b), 
ATLAS pj spectrum data on {Z/y*) production to (bare) e + e~ pairs, the circular dots are the data, 
the blue(green) lines are HERWIRI1.031(HERWIG6.510). In both (a) and (b) the blue(green) squares 
are MC@NLO/HERWIRI1.031(HERWIG6.510(PTRMS = 2.2GeV)). In (b), the green triangles are 
MC@NLO/HERWIG6.510(PTRMS =0). These are otherwise untuned theoretical results. 

analysis in Ref. [9] of the FNAL data on the single Z/y* production in pp collisions at 1.96 TeV. 

In Fig. 1 1 of the second paper in Ref. [9], we showed that, when the intrinsic rms pj parameter 
PTRMS is set to in HERWIG6.5, the simulations for MC@NLO/HERWIG6.510 give a good fit 
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to the CDF rapidity distribution data [28] therein but they do not give a satisfactory fit to the DO 
p T distribution data [29] therein whereas the results for MC@NLO/HERWIRI 1.031 give good fits 
to both sets of data with the PTRMS = 0. Here PTRMS corresponds to an intrinsic Gaussian 
distribution in pj. The authors of HERWIG [27] have emphasized that to get good fits to both 
sets of data, one may set PTRMS = 2 GeV. Thus, in analyzing the new LHC data, we have set 
PTRMS = 2.2GeV in our HERWIG6.510 simulations while we continue to set PTRMS = in our 
HERWIRI simulations. 

Turning now with this perspective to the results in Fig. ||, we see confirmation of the finding 
of the HERWIG authors. To get a good fit to both the CMS rapidity data and the ATLAS p T data, 
one needs to set PTRMS = 2GeV [30] in the MC@NLO/HERWIG6510 simulations. We again 
see that at LHC one gets a good fit to the data for both the rapidity and the pj spectra in the 

MC@NLO/HERWIRI1.031 simulations with PTRMS = 0. In quantitative terms, the % 2 /&.o.f. for 
the rapidity data and p T data are (.72,.72)((.70,1.37)) for the MC@NLO/HERWIRI1.031(MC@NL 
O/HERWIG6510(PTRMS=2.2GeV)) simulations. FortheMC@NLO/HERWIG6510(PTRMS=0) 
simulations the corresponding results are (.70,2.23). 

The usual DGLAP-CS kernels require the introduction of a hard intrinsic Gaussian spread 
in p T inside the proton to reproduce the LHC data on the p T distribution of the Z/y* in the pp 
collisions whereas the IR-improved kernels give a better fit to the data without the introduction of 
such. This hard PTRMS is entirely ad hoc; it is in contradiction with the results of all successful 
models of the proton wave-function [31], wherein the scale of it is < .4GeV. More importantly, 
it contradicts the known experimental observation of precocious Bjorken scaling [32,33], where 

the SLAC-MIT experiments show that Bjorken scaling occurs already at Q 2 = 1+ GeV 2 for Q 2 = 

—q 2 with q the 4-momentum transfer from the electron to the proton in the famous deep inelastic 
electron-proton scattering process whereas, if the proton constituents really had a Gaussian intrinsic 
pj distribution with PTRMS = 2GeV, these observations would not be possible. What can now 
argue is that the ad hoc PTRMS = 2.2GeV value is a phenomenological representation of the more 
fundamental dynamics realized by the IR-improved DGLAP-CS theory. Is it possible to tell the 



lould be able to distinguish the two. 



difference between the two representations of the data in Fig. 

Physically, one expects that more detailed observations sL . _ 

Indeed, in Fig. | we show the MC@NLO/HERWIRI 1.031 (blue squares) and MC @ NLO/HER 
WIG65 10(PTRMS=2.2GeV) (green squares) predictions for the Z/y* mass spectrum when the 
decay lepton pairs are required to satisfy the LHC type requirement that their transverse momenta 

{p T ,Pj} exceed 20 GeV. As the peaks differ by 2.2%, the high precision data such as the LHC 
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Figure 3: Normalized vector boson mass spectrum at the LHC for pr (lepton) > 20 GeV. 
ATLAS and CMS experiments will have (each already has over 5 x 10 6 lepton pairs) will allow 
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one to distinguish between the two sets of theoretical predictions. Other such detailed observations 
may also reveal the differences between the two representations of parton shower physics and we 
will pursue these elsewhere [12]. In closing, two of us (A.M. and B.F.L.W.) thank Prof. Ignatios 
Antoniadis for the support and kind hospitality of the CERN TH Unit while part of this work was 
completed. 
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